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With the discovery of massive neutron stars such as PSR J1614-2230, the question has arisen 
whether exotic matter such as hyperons can exist in the neutron star core. We examine the conditions 
under which hyperons can exist in massive neutron stars. We consistently investigate the vector 
meson-hyperon coupling, going from SU(6) quark model to a broader SU(3) symmetry. We propose 
that the maximum neutron star mass decreases linearly with the strangeness content f s of the 
CN| . neutron star core as M ma x{fs) = M ma x(0) — O.6M0(/ S /O.l), which seems to be independent of the 

underlying nuclear equation of state and the vector baryon-meson coupling scheme. Thus, pulsar 
mass measurements can be used to constrain the hyperon fraction in neutron stars. 
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I. INTRODUCTION 

The discovery of the massive neutron star PSR J1614-2230 has raised new challenges for theories of dense matter 
beyond nuclear saturation density. Shapiro delay measurements from radio timing observations of the binary 
millisecond pulsar indicate a large mass of 1.97±O.O4M of the neutron star The core of a neutron star harbors 
(i a dense matter environment, which could be the site for strangeness containing matter, such as hyperons. Though 
nuclear interactions in the saturation regime are well understood, one has to utilize neutron star observations to find 
+3 , clues about the physics of cold and dense matter beyond several times saturation density. Any theory of ultra dense 
^ ■ matter has to explain the recently observed large neutron star mass. According to existing models of dense matter, 
i '_ the presence of hyperons leads to a considerable softening of the equation of state (EoS), resulting in a corresponding 
reduction of the maximum mass of the neutron star [2j-l5[- Then the existing theories involving hyperons are in 
' conflict with the large pulsar masses On including hyperons, most relativistic models obtain maximum neutron 
star masses in the range 1.4 — 1.8M 0413] ■ However, in exceptional cases, neutron stars with maximum masses 
1 larger than 2Mq have been obtained, either by pushing the threshold for appearance of hyperons to higher densities, 
. or due to strong hyperon vector repulsion (l#{2l1 |. Taurines et al. [22| achieved large neutron star masses including 
■ hyperons by considering a model with density dependent coupling constants, which were varied nonlincarly with the 
scalar held. Recently, Bednarek et al. [lp also achieved a stiffening of the EoS by using a non-linear relativistic mean 
field (RMF) model with quartic terms involving the hidden strangeness vector meson. In addition to the inclusion 
of such a meson into a density dependent RMF model, Lastowiecki et al. [24[ assumed a quark matter core in order 
to obtain massive stars. Bonanno and Sedrakian [25j also modeled massive neutron stars with a hyperon and quark 
core using a fairly stiff EoS and vector repulsion among quarks. In several studies, the maximum neutron star masses 
obtained when including hyperons were not very different from those containing nucleons only [HI UtI [l9j . In more 
' sophisticated models such as the Brueckner-Hartree-Fock model, the maximum neutron star masses were generally 
found to be lower than 1.6M Q which is in contradiction with observed pulsar masses [3j-[ol f26l - [30j . 
From the studies cited in the previous paragraph, it seems that the possible presence of hyperons in massive neutron 
stars is in many cases reconciled by incorporating large vector repulsion in an ad hoc way. In contrast, we investigate 
the role of vector repulsion starting from symmetry arguments. Assuming SU(3) symmetry, we perform a controlled 
parameter study and constrain the parameters using the observed mass of PSR J1614-2230. This procedure is in 
line with modern microscopic models for realistic baryon-baryon potentials such as the Bonn potentials (3l| and 
the Nijmegen potentials [331 which adopt SU(3) symmetry to describe the baryon interactions for the baryon octet. 
For our investigations, we employ a RMF model, in which the parameters are calibrated around nuclear saturation 
density 0, Hj]. However, the extrapolation of such properties to supranuclear densities presents uncertainties. In 



'Electronic address: s.wcisscnborn@thphys.uni-hcidclbcrg.de 
^Electronic address: d.chatterjee@thphys.uni-hcidclbcrg.de 
t Electronic address: schaffncr@thphys.uni-heidclberg.de 



2 



a previous paper [341 ] . we investigated how the uncertainty in nuclear saturation properties, such as effective nucleon 
mass or nuclear compressibility, or hypernuclear properties, such as potential depths of hyperons in nuclear matter, 
could influence our conclusions about the presence of hyperons in the core of massive neutron stars. In this work, 
we question the fundamental assumption of SU(6) symmetry, which relates the hyperon couplings to the nuclear 
couplings. 

This paper is organized in the following way. In Sec. HH we describe the model to calculate the EoS. The parameters 
of the model are listed in Sec. IIII1 The results of our calculations are discussed in Sec. IIV1 and a summary of our 
conclusions is given in Sec. [Vj 



II. THEORETICAL MODEL 



As elaborated in our previous paper [34|, a RMF theoretical model is adopted to describe neutron star matter 
subject to chemical equilibrium and charge neutrality The octet baryons, electrons and muons are considered as 
constituents of the core. In this model [331 ] . baryon-baryon interaction is mediated by the exchange of scalar (a), 
vector (lo) and isovector (p) mesons. The hypcron-hyperon interaction is incorporated through additional strange 
scalar (a*) and vector ((/>) mesons. The stiffest possible EoS within the model is obtained on inclusion of the strange 
vector meson <fi and by omitting the a* meson, and is referred to as the "model cro;p0" [34|. The vector meson- hyperon 
couplings in this model are related to those of nucleons through the symmetry of the SU(6) quark model. In this 
paper a consistent investigation for the vector meson coupling is achieved, going from the SU(6) symmetry to the 
more general SU(3) symmetry, in order to provide a general analysis for the role of repulsive hyperon interactions on 
the properties of neutron stars. 



A. Determination of vector meson-hyperon couplings 



1. Flavor SU(3) and the baryon-meson interaction Lagrangian 



SU(3) in flavor space can be regarded as a symmetry group of strong interaction (restricting it to three quark 
flavors, up, down and strange). For neutron stars, we only consider the [8] representation for the baryons, namely the 
J p = i octet. The product of the baryon and meson representations reads: 

[8] <g> [8] = [1] [8} A © [8] s © [10] © [27]. (1) 

With the help of matrix representations for the baryon octet (B) and meson nonet (singlet state Mi and octet states 
Ms), SU(3) invariant expressions can be constructed (32|. The interaction Lagrangian for the whole meson nonet and 
the baryons can be written as a sum of terms, one coming from the coupling of the meson singlet to the baryon octet 
(S), and other two terms from the interaction of the meson octet and the baryons - one being antisymmetric (F) and 
the other being symmetric (D) [35L l36l] : 

C mt = -g s V2 [aTr ([B, M 8 ]B) + (1 - a)Tr ({B, M 8 }B)] - g x J^Tr^B^^) . (2) 

Here, g% and g\ denote the meson octet and singlet coupling constants. The F/(F + D) ratio a, which lies in the 
range of < a < 1, is a weight factor for the contributions of the symmetric D and the antisymmetric F couplings 
relative to each other. 

The assumption of SU(3) symmetry implies that all possible combinations of couplings for each type of meson with 
all possible baryons can be described with only four parameters: the singlet coupling constant g-y, the octet coupling 
constant g$, the F/(F + D) ratio a [HI, H3] and a "mixing angle" 9 relating the physical isoscalar mesons to their 
pure octet and singlet counterparts. The SU(6) quark model is obtained as the combination of flavor SU (3) with 
spin SU(2) (and is thus a special case of flavor SU(3 j). 

We restrict our attention to vector coupling terms containing only baryons of the same species [38| and denote the 
F/(F + D) ratio corresponding to vector meson interaction by ay. The requirement of spin independence for the 
couplings of the identically flavored £ and A baryons: 



gsui = gAu, gs<t> = gA$ 



(3) 
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leads to a fixed value of the Fj (F + D) ratio ay: 

a v = l . (4) 

This means, we have a pure i^-type coupling in the interaction Lagrangian ([2]). Further, the proposition that the 
nucleon does not couple to the 0- meson which is a pure ss-state (i.e. gN<f, = 0), in the case of "ideal mixing" for the 
ui and (f> meson (39j . 

tanOy = —j= , (5) 
v2 

leads to the relative coupling strengths [13 

gNu ■ gAu ■ g^u : gsw =3:2:2:1 (6) 

and 

3A0 : : = 1 : 1 : 2 (7) 
where equations © and (J7J are related through [38| : 

.9A^ = --^ffH0 • (8) 

This leaves only one degree of freedom, say gNu which is fitted to saturation properties of nuclear matter Q. We fit 
the a meson hyperon couplings to the potential depths of hyperons in nuclear matter (see [34[ for discussions) . 



2. Beyond SU(6) 

In order to explain the large observed mass of PSR J1614-2230 we might either introduce more particles and 
parameters into our RMF model, or reconsider the necessity of fixing them to the SU (6) values. 

From the quadratic mass formula for the mesons one obtains 6y « 40° which is quite close to the ideal mixing 
angle 9 ~ 35.3° [39| . Thus, we may retain the condition of ideal mixing for the vector mesons. The parameters ay, 
gi and gs however, we may consider as being free. Since we later fix gNu to saturation properties of nuclear matter, 
we combine the singlet and octet coupling constants to a single parameter z which we define as 

z := f (9) 
.91 

and keep the parameters ay, z and g N ^„ If we set z to its SU(6) value z = 1/ s/Ti and use ideal mixing while varying 
ay, which by definition lies within the range < ay < 1, it gives: 

gAu _ 2ay + 4 gzu _ 8 - 2ay g E u _ 5 - 2ay 
gNu 4ay + 5 gNu 4ay + 5 gNu 4ay + 5 

g M = v __4ay-4 
gNu 4ay + 5 

gAj, _ /j) 2ay - 5 gsj, _ 2a v + 1 g E4> _ ^ 2a v + 4 



gNu 4ay + 5 gNu 4ay + 5 gNui ^Q-v + 5 

If we set ay to its SU(6) value ay = 1 and use ideal mixing while keeping z as a free parameter instead we arrive at: 

gAui _ \/% _ gsu gsu _ V% - \/3z 

gNu y/2 + Viz gNu gNu s/2 + s/3z 

gN4> _ V&z - 1 g A(j> _ -l _ gs(j> g~<t> _ l + VEz 

gNu \/2 + Viz gNu \/2 + \/3z gNu gNu y/2 + y/3z 
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If we require the interaction due to u exchange to be repulsive for all baryons, we want no changes of sign in the ui 
couplings, i.e. especially gsu/dNu > and therefore < z < 2/V6- The couplings of the p meson to the baryons 
may be related to g^u by relations analogous to (|T0|) . ([TT|) involving ay or z. However, since the p couplings control 
the asymmetry energy and its density dependence L, g^ p may alternatively be fitted directly to the asymmetry 
energy coefficient at saturation (following Ref. [loj]). For the hyperons, the p coupling strengths are then related to 
gNp according to isospin. This method we adopt for our calculations and thereby we avoid the problem of having 
gNp = at z = 0, ay = 1 which would yield unphysically low values for the symmetry energy coefficient. For the 
two ways to fix the p couplings, the difference in the obtained maximum neutron star masses is rather negligible 
(AM max < O.O5M ). 

From the relations (|10|) and it is clear that for values of ay or z not coinciding with SU(6) the <f> meson couples to 
the nucleon which is supported by the large strange quark condensate in the nucleon found in lattice gauge simulations 




III. PARAMETERS OF THE MODEL 
A. Nucleon-Meson coupling constants 



The nucleon- meson coupling constants g„N, g^N, 5pv, b and c are determined from the saturation properties of 
nuclear matter [40| : binding energy = —16.3 MeV, baryon density no = 0.153 fm -3 , and asymmetry energy coefficient 
Oasy = 32.5 MeV. The incompressibility K and effective nucleon mass m* N /tun are varied according to the parameter 
set in consideration (e.g. GM1 @, NL3 H, TM1 Q). 



B. Hyperon-Meson coupling constants 



The strange and non-strange vector meson-hyperon couplings have already been described in detail in the previous 
section. The non-strange scalar meson-hyperon couplings are fitted to the potential depths for the hyperons in nuclear 
matter. The following values l/[ N) = -30 MeV, U^ N) = +30 MeV, ui N} = -28 MeV Jiave been adopted from the 



hypernuclear experimental data (see e.g. [4514481 ] , a summary of usual choices is found in (49|). However, the particular 
choice of hyperon potentials does not have crucial consequences regarding the maximum mass of neutron stars as has 
been discussed in [34j . In fact, a larger S potential is the more conservative choice as it gives a slightly smaller 
maximum mass. 

When fixing the parameter g^^ we remember that by means of non- vanishing g^^ the (j) meson also contributes 
to the saturation properties of nuclear matter. In particular, in all thermodynamic quantities, e.g. in the nucleon 
chemical potential, we have to replace the term 

gNujti — > Siv^w + 9N4>4> (12) 

where the tilde ( ~ ) denotes the coupling and corresponding field at nuclear saturation for the case of gN<j, = 0. RMF 
models without vector meson self- interaction (e.g. GM1, NL3) yield as equations of motion for the w and cj>: 

g Nu) u = — T p N , g Nrj) (p = — T p N . (13) 
mi m% 

This allows to rewrite the substitution prescription (|12[) as 

9 2 N^9%Jl + § l 3) ( 14 ) 



9nu 



where the ratio gN<t>/ gNu> appearing inside the term in parentheses is a function of ay or z as it fulfills equation (|10[) 
or (|TT]) respectively. We fix g^^ to the saturation properties of nuclear matter (following Ref. [ioj ) , because then it is 
clear that the EoS of pure nuclear matter is independent of ay and z. This is desired, because with or without <f> we 
wish the properties of nuclear matter to be independent of any variations in the hyperon coupling strength relative to 
the nucleon one. In the TM1 model there is an additional d(w A ') 4 /4 term in the Lagrangian that prevents the mean 
meson fields from droping out of the equation when substituting as in (fT2|) . Thus, it becomes necessary to solve the 
nonlinear equation of motion for the new and old u at nuclear saturation density. To fully resolve the problem of 
nonlinear vector meson self-interactions we would need to adapt g^u, as well as d to the presence of (f> mesons, which 
requires an additional input from experiments. 
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FIG. 1: Relative vector meson coupling constants as functions of the gs/gi ratio z for fixed ay = 1. The value z = 1/^6 
corresponds to the SU(6) case. 



IV. RESULTS AND DISCUSSIONS 

A. Varying the gs/gi Ratio z 

We probe the effects of the gs/gi ratio z on the stiffness of the hadronic EoS. We plot the coupling constants as 
functions of z in Fig. [TJ for the explicit formulae given in equation (jll[) . For z = 0, all coupling constants gsu are 
the same and similarly are all coupling constants gB<f> equal. This is due to the fact that z = corresponds to g$ = 0, 
which results in the equality of the corresponding baryon-meson couplings as the baryons couple only to the flavor 
singlet state. With increasing z, i.e. with increasing contribution from the coupling to the octet g%, the resulting 
couplings all become smaller except for gs<t>- At z = 1/V& ~ 0.4082 the SU(6) case is reached where the </> does 
not couple to the nucleon. Thereafter, for z > l/-\/6 the coupling constant g^^ changes its sign so that it is not a 
repulsive but now an attractive interaction. Note, that the uj coupling constants for A and £ hyperons are equal for 
all values of z. As anticipated in Sec. Ill A 21 we restrict z to the interval z G [0 : 2/v / 6] . 

The EoS for z = 0.1, 0.2, 0.8 are plotted in Fig. [5J At first glance it becomes clear that the EoS stiffens with 
decreasing z. This can be explained with the help of Fig. Q] where we had plotted the z dependence of the vector 
meson coupling constants. We noticed that with increasing z, all couplings except gs<f> decrease. Since the S hyperons 
only play a subordinate role compared to the neutrons, the increase of gs<j> does not prevent that part of the overall 
interaction between the baryons, which is mediated by the vector mesons, to become less repulsive. Therefore, the 
EoS must soften with increasing z. Together with further decreasing coupling strengths between the vector mesons 
and the other baryons (except again g~<j>) the EoS becomes even softer until z = 1j\f% is reached. 

In Fig. [3] we plot the mass-radius relations for the various EoS we just discussed in the context of Fig. [21 As 
expected from the influence of z on the stiffness of the EoS, the lowest maximum mass is obtained for z = 2/\/6 sa 
0.8165, or in the case of Fig. [3] at z = 0.8, namely M = 1.49A/q. The maximum mass grows up to the value 
M = 2.36-Mq for z = 0. We notice that the maximum mass of a neutron star in our RMF model reacts rather 
strongly to the variation of z: over the whole z range the change in the maximum mass is AM=0.87Af Q . 

After varying z in rather big steps, we now plot in Fig. [4] the maximum neutron star mass as a continuous function 
of z for the models GM1, NL3 and TM1 and for nuclear matter as well as for baryonic matter. As already analyzed in 
the discussion of Fig. [TJ we see that the branch for nucleonic matter is insensitive to the changes in z for the NL3 and 
GM1 models. For TM1, the quartic self interaction term in the Lagrangian spoils this property since we had to keep 
the corresponding coupling constant d from the SU(6) value of z fixed for the whole z range: the maximum masses 
for pure nucleonic stars therefore depend on the actual z value, showing a minimum for the SU(6) case z w 0.4082, 
namely M = 2.18Af©. 

Considering the neutron stars containing hyperons, we see in Fig. 0] that the maximum masses depend on z as 
already observed in Fig. [3] for the largest z values the maximum masses are the smallest and continually grow with 
decreasing z. It is interesting to see that towards z — the maximum masses of these stars seem to approach the 
maximum masses of the corresponding pure nucleonic stars for all parameter sets studied. A look at the particle 
number fractions for the GM1 parameter set at z = 0.8 (which we plot in Fig. E^a)) shows that the first hyperons 
to appear in the hadronic matter are the 5~ and the A at total baryon number densities of nf, rj 0.28 /m~ 3 and 
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FIG. 2: EoS for different gs/gi ratios z within a nonlinear a — uj model with additional <j> meson and the full baryon octet for 
GM1 parameter set. The EoS get stiffer with decreasing z. 

rif, ~ 0.29 ,fm~ 3 respectively, while the S° appears much later at rib ~ 0.76 fm~ 3 . On increasing z to its SU(6) value 
(Fig. 15(b)), the A hypcron appears first at rib ~ 0.36 fm~ 3 , followed by S~ at rib ~ 0.4 fm~ 3 and S° at rib ~ 0.89 
fm~ 3 . At z = (Fig. 15(c)), the threshold of appearance of hyperons is pushed to even higher densities: rib ~ 0.73 
fm~ 3 for A, rib ~ 0.74 fm~ 3 for S~ and rib ~ 1-38 fm~ 3 for 3°. Thus, for z = the neutron stars consist mainly of 
nuclear matter which is why the maximum masses are so close to those of the pure nuclconic stars. For the case of 
the parameter sets NL3 and TM1, the particle fractions are qualitatively the same as in the GM1 case. In the case 
of N L3 p arameterization, a well-known instability occurs at high densities when the effective nucleon mass becomes 
zero [331 ] . The critical density for the appearance of the instability depends on the value of the hyperon coupling 
constants. However for the present investigation, this instability plays no role as it appears beyond the maximum 
densities reached in the neutron star interior. 



1. Combining m* N and z variations 

The impact of z on the maximum mass of neutron stars is as comparably large as the influence of the effective 
nucleon mass at saturation m* N as investigated in our previous study [34| . We therefore combine both parameters in 
a single plot, Fig. [6j where we show the maximum neutron star mass as a function of m^/mjv for different z values. 
The incomprcssibility in this case is fixed to K = 240 McV, but the exact value is irrelevant as shown in our previous 
paper [34| . We see in Fig. [6] that the effective mass has basically the same effect for all z values and z the same 
effect for all effective masses: for fixed z, the maximum masses decrease drastically for increasing effective mass. For 
low z values, where the EoS is stiffer than for higher z values, the dependence of the maximum mass on the effective 
mass is slightly larger: the difference along the whole range of m* N /mM is ps 0.6M© for z = 0.8 while for z = it is 
w 0.9Mq. The influence of z on the maximum masses is slightly more pronounced than that of the effective masses: 
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FIG. 3: Mass-radius relations for the EoS displayed in Fig. 
baryons couple to the vector mesons with equal strengths. 



[2] The maximum mass is obtained for the case z = 0, where all 
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FIG. 4: Maximum masses as functions of the gs/gi ratio z for NL3, GM1 and TM1 parameter sets. For each parameterization 
the case of pure nucleonic matter is also displayed. 




FIG. 5: Particle fractions for the GM1 parameter set for 3 different values of z: z = 0.8 (a), z = 0.408 (b) corresponding to 
the SU(6) case, and z = (c). The threshold for the appearance of the hyperons A, E~ and H° is pushed to higher densities 
with decreasing z. 

the difference in the maximum mass between z = and z = 0.8 is ss 0.65M Q for m* N /m^ = 0.8 and about IMq for 
m* N /mN — 0.55. For comparison, we also plot the maximum masses of purely nucleonic neutron stars, and also mark 
in the figure the points corresponding to the SU(6) case, for several other RMF sets fitted to properties of nuclei e.g., 
TM1 , NL3 or NL-SH [H, 0, [50|-[52j . In Fig. 3] the masses of baryonic stars were found to approach the limit of purely 
nucleonic stars for decreasing values of z. With Fig. [5] we can now visualize that the maximum masses in these two 
cases differ slightly (A M < O.OIMq). 

We close this section by concluding that a maximum mass of at least 1.97 ± O.O4M0 requires very small values of z 
for large effective masses (z < 0.1 at m* N /rriN = 0.8), z values around SU(6) for effective masses close to m* N /rriN w 0.7 
and very low effective masses for large z values (m^/mjv < 0.58 for z = 0.7). We note, that z values close to the 
maximum of z = 2/V6 are now allowed configurations within the plot range: the maximally allowed z value for the 
investigated model is z < 0.77 at m* N /mM = 0.55. 

B. Varying the F/(F+D) ratio ay 

After systematically investigating a wider z range we repeat the formalism for the ay-ratio in the present section. 
The F/(F + D) ratio is by definition restricted to the interval ay € [0; 1], where the lower bound corresponds to a 
pure D-type coupling and the upper limit (i.e. the SU(6) value) corresponds to a pure F-type coupling. Analogous to 
the case studied above, we adopt ideal mixing as well as a g%/g\ ratio fixed to its SU(6) value z = l/v6 and we allow 
for the 4> meson to couple to the nucleon. We plot the coupling strengths in Fig. [7] where we vary ay between and 
1. Note that the coupling constants <? u a for A hyperons do not change considerably and that all vector couplings 
remain repulsive (do not change their sign). 

The ratio of the F- and D-type couplings can be continuously varied between the two extremes of a pure F- and a 
pure D-type coupling. We see that from ay = 1 down to ay = all couplings become stronger except for gs^/gNui- 
Since the £ hyperons have but very little influence on the EoS up to neutron star densities, we can expect that the 
EoS become stiffer with decreasing ay. This is exactly what we find in Fig. [8] where we plot the EoS for several 
values of ay using "model au>p4>" and GM1 parameter set. 

The stiffness depends monotonously on ay and we get the softest EoS for the SU(6) case ay = 1, i.e. a pure 
F-type coupling, while the stiffest EoS is obtained for ay = which corresponds to the pure D-type coupling of 
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FIG. 6: Maximum masses of hyperonic neutron stars as functions of effective nucleon mass m* N for different values of the gs/gi 
ratio z. For comparison, a line for nucleonic stars and points to mark RMF sets (e.g. TM1, NL3) corresponding to the SU(6) 
case are also given. 
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F/(F+D) ratio a v 



FIG. 7: Vector meson coupling constants as functions of the F/(F + D) ratio ay. The gs/gi ratio is fixed to its SU(6) value 
z = l/\/6 and ideal mixing is assumed. The SU(6) case is given by av = 1. 



the baryon and meson multiplcts. The EoS for ay < 0.2 appear to be indistinguishable at neutron star densities. 
This is evident in Fig. [9j where we plot the mass-radius relations corresponding to the EoS from Fig. [8j for the 
values ay = 0.0 — 0.2 the maximum masses but also the radii of the corresponding stars coincide (M max = 2.36M Q , 
R= 11.8 km). We note, that this value of the maximum mass is also obtained for the purely nucleonic case (compare 
e.g. Fig. [4j. Thus, for the "model auipcf)" the nuclear matter limit is reached below ay < 0.2 in the case of GM1, 
GM3 and NL3 parameter sets, while for the very stiff PL-Z EoS (having effective mass m* N /mN — 0.55) [5l| pure 
nucleonic stars are already obtained for ay < 0.3. 

In this way, for hyperonic stars the limit of nucleonic stars is continuously approached for decreasing values of the 
gs/gi ratio z, or for decreasing values of the F/(F + D) ratio ay away from the SU(6) value, respectively. 

To generalize our findings, we plot in Fig. [10] the maximum masses of neutron stars as a function of 
strangeness fraction f s (the number of strange quarks divided by the total number of quarks) for four different EoS: 
for m^/mN = 0.55, 0.75 and for the NL3 (m* N /rnN = 0.6) and GM1 (m* N /rnN = 0.7) parameter sets, by varying 
ay (solid lines) and z (dashed lines). It is evident from the figure that on decreasing ay or z, the strangeness 
fraction in the core decreases, and there is a corresponding increase in the maximum mass of the star. At zero 
strangeness fraction the nucleonic limit is reached, and this corresponds to the highest value of the maximum mass. 
For m* N /mN ~ 0.7, the relation between the maximum mass of the star and its strangeness fraction can be fitted 
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900 




FIG. 8: EoS for "model aujp4>" in GM1 parameterization for different values of oiv ■ z is fixed to its SU(6) value z = 1/y/Q and 
ideal mixing is assumed. The EoS become stiffer with decreasing ay. 



linearly according to the formula: 

M max M max (f s = 0) / f s \ 

~M^ = ^ C lcuJ' (15) 

where c ~ O.GMq. This has interesting consequences when we use these results to predict the maximally allowed 



strangeness fraction in maximum mass neutron stars. In an associated study [531 ] . we applied the results from 
measurements of sub-threshold kaon production in heavy-ion collisions (HIC) to study the implications on neutron 
star properties. It was found that the heavy-ion data and causality imply a firm upper limit on maximum mass 
of compact stars of 3 solar masses. Substituting this value into the derived formula above, we can show that the 
strangeness fraction in a 2 M Q star cannot be more than: 



M HIC - M obs 

fmax max max 

Js 



6M Q 
6Mm 



0.17 . (16) 



We also point out that for the EoS compatible with the observed mass limit, the hyperon fraction in a canonical 
IAMq star is zero. Only for small effective nucleon masses and values of z above the SU(6) value can a hyperon 
fraction of less than 0.5% be reached. 

We have now found that we can reach the nuclear matter limit starting from a RMF model including hyperons and 
continuously changing the model parameters, instead of making a discrete "on/off" decision about whether to include 
or exclude hyperons. Instead, we should take the position of saying that one can only exclude a RMF parameter set 
as soon as the corresponding maximum mass for nucleonic stars is incompatible with observations. As long as the 
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FIG. 10: Maximum masses of neutron stars as functions of strangeness fraction f 3 in the neutron star core for four different EoS. 
The maximum mass decreases linearly with the strangeness fraction approximately as Mm ax (f a ) = M max {Q) — 0.6M©(/ S /0.1). 
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observational mass limit is below the nucleonic mass limit of the model, it is possible to have hyperons in the core of 
maximum mass neutron stars. 



V. SUMMARY 



We investigated the conditions for the existence of hyperons in massive neutron stars. We went beyond the 
spin- flavor SU(6) quark model for determining the vector meson couplings to the more general relations from flavor 
SU(3), where we varied the g%/gi ratio z and the F/(F+D) ratio ay while assuming ideal mixing. 

We fixed either z or ay to their SU(6) values and allowed for a non-vanishing N — <f) coupling. We found within 
"model cruip<f> v , that decreasing z below its SU(6) value leads to stiffcr EoS and larger maximum masses. The most 
massive stars were obtained for z = 0, where only the meson singlet couples to the baryon octet, meaning that all 
vector couplings take the same value. We have seen that at z = the masses of the stars come very close to the case 
of purely nucleonic neutron stars. 

Keeping z fixed at its SU(6) value and varying ay has even a greater impact on the maximum masses of neutron 
stars: the lowest value is given for the SU(6) case ay = 1. Towards ay = the maximum masses increase 
monotonically. In both cases, i.e. on lowering z and on lowering ay, we found that the combined repulsive 
interactions mediated by the u and <f> mesons cause the hyperons to appear at successively higher densities so that 
hyperons appear only in a small core (z = 0, ay = 1) of the maximum mass star or not at all (z = 1/V6, ay = 0). 
Thus, the EoS becomes stiffer and the strangeness fraction of neutron stars decreases with decreasing values of ay 
or z until the limit is reached where the whole sequence of neutron stars contains only pure nucleonic stars. There is 
a smooth transition between strange hadronic and nucleonic stars as well as between their corresponding maximum 
masses where formerly there was only an explicit inclusion/exclusion of hyperons in the construction of the model. 
This finding holds for all investigated parameter sets, but it should be stressed that the inclusion of the 4> is vital for 
reaching the nuclear limit. Wc find that the maximum mass decreases linearly with the strangeness fraction f s of the 
neutron star core as M max (f s ) = M max {Q) — 0.6Mq(/ s /0.1), independent of the chosen nuclear EoS and the adopted 
vector baryon-meson coupling scheme. For instance, given the mass of 3Mq for a purely nucleonic neutron star, the 
strangeness fraction for a neutron star with an observed mass of 1.97M Q would then be below f s j§ 17%. 
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